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Abstract: Starting from Sklyanin's separation of variables for the sl^ Yangian model, we derive 
the separation of variables for the quantum 5^3 Gaudin model. We use the resulting new variables 
for rewriting the si^. Knizhnik-Zamolodchikov equations, and comparing them with certain null- 
vector equations in conformal field theories with Wa-algebra symmetry. The two sets of equations 
are remarkably similar, but become identical only in the critical level limit. This is in contrast to the 
Knizhnik-Zamolodchikov equations, which are known to be equivalent to Belavin-Polyakov- 
Zamolodchikov equations for all values of the level. 



Contents 



|l[ Introduction and conjecture 

^ Gaudin Hamiltonians in conformal field tlieory 

[2.1| Knizhnik-Zamolodchikov equations 
O Ward identities for the cubic field 



^ Review of the si2 case 

p.l| Separation of variables for the s£2 Gaudin model 

p.2| The s^2 Knizhnik-Zamolodchikov equations in Sklyanin variables 

^3] Comparison with Virasoro null-vector equations 



Tlie case JT] 
[4.1| Separation of variables for the 5^3 Gaudin model |ll] 
[4.2| The s^s Knizhnik-Zamolodchikov equations in Sklyanin variables 14 



3] W3 null-vector equations 16 
]4| Comparing 5^3 Knizhnik-Zamolodchikov equations with W3 null-vector equations 20 



^ Conclusion 22 

^ A few teclinical results 23 

|Al| Helpful identities ^ 



A.2| A characterization of ^2 (y) [24 



A.3| Study of the differential operator Di 24 



1. Introduction and conjecture 

Many interesting models of two-dimensional conformal field theories are based on affine Lie alge- 
bras s£n- and their cosets, starting with Wess-Zumino-Witten models. To solve such theories is an 
interesting challenge, whose difficulty depends more from the choice of the underlying Lie algebra 
s£n, than from the particular coset or real form chosen. 

For example, the 5^2 family includes string theory in AdSs and in the 5'L(2,]R)/f7(l) 2d 
black hole, as well as the model; the simplest non-rational nontrivial model of the family 
is however Liouville theory, also known as conformal 5^2 Toda theory. In several of the other 
theories in the s£2 family, it turns out that arbitrary correlation functions have a simple relation 
to certain Liouville theory correlation functions This relation entails a relation between the 

Knizhnik-Zamolodchikov equations which follow from 5^2 symmetry, and the Belavin-Polyakov- 
Zamolodchikov equations which follow from the conformal symmetry of Liouville theory [^. The 
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relation to Liouville theory is helpful in solving certain models in the sH.2 family, by disentangling 
the particular details of a model from its general s£2-based properties. For example, the H^- 
Liouville relation was very helpful in solving the model on a disc [Q]. Moreover, playing 
with the Liouville side of the relation leads to the discovery of new conformal field theories which 
generalize the model [|5|], and which can be considered as members of an extended 5^2 family. 

The intuitive reason why such a relation exists is that 5^2 representations are parametrized by 
just one number, their spin. So it is not very surprising that the dynamics of say the model, 
a theory of three interacting bosons, are in some sense effectively one-dimensional. Applied to 
a theory with an s£n>2 symmetry algebra, which may involve as many as — 1 bosons, this 
reasoning suggests that it could be related to a theory of only — 1 bosons. Such a theory is 
present in the si^ family: namely, conformal si^ Toda theory, which can be described by the 
Lagrangian L = {d<j), dcp) + eK^i^'l') where the field (j){z, z) and the simple roots live in 

the — 1-dimensional root space of sl^. (See for example ^ for details.) It is therefore natural 
to investigate whether correlation functions of that theory have a simple relation to correlation 
functions of other models in the family. Such a relation would be a welcome simplification: for 
instance, in the 5^3 family, we would trade 8 bosons of the 5'L(3, M) WZW model for the 2 bosons 
of slz conformal Toda theory. 

The investigation of the sIn>2 families is motivated both from the appearance of groups of 
rank higher than one in many interesting string theory backgrounds, and from the observation that 
theories in the s£n>2 families are qualitatively more difficult, and more generic, than theories in the 
s^2 family. This is due to features like: infinite fusion multiphcities, correlation functions involving 
degenerate fields without obeying nontrivial differential equations, and structures constants which 
can probably not be written in terms of known special functions [^]. These are serious obstacles 
in the way of solving such theories. Nevertheless, we do know a strong explicit constraint on the 
correlation functions of all models which have the full sIn symmetry: they obey KZ equations. 
The aim of the present article is therefore to determine whether the sis KZ equations are related to 
some null- vector equations in conformal 5^3 Toda theory, which follow from its symmetry algebra 

m. 

In analogy with the 5^2 case, we will look for a relation based on Sklyanin's separation of 
variables [0]. As the KZ equations are closely related to the Gaudin Hamiltonians, we will use 
Sklyanin's separation of variables for the quantum sis Gaudin model. Before using it, we will 
actually have to work it out, as this has apparently not been fully done in the existing literature. A 
rather close starting point is available though: the separation of variables for the sis Yangian model 

Let us now sketch the correlation functions we are interested in and the relation we are aiming 
at. Consider a theory with an 5^3 symmetry algebra at level k. We are interested in correlation 
functions of generic 5^3 affine primary fields ^^{x\z), where the spin j labels 5^3 representations, 
the variable x is a generic isospin coordinate (a triplet of complex numbers), and z is a coordinate 
on the complex plane where the field lives. We denote an n-point function of such fields as 

nn=(j{^'^{x,\z,)^ . (1.1) 

We will seek to relate such correlation functions to fairly particular correlation functions in a theory 
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with a Ws symmetry algebra at parameter 6 = (A; — 3)~2 , which involve not only n generic W3- 
primary fields Va^{zi) corresponding to ^^^{xi\zi), but also 3n — 6 degenerate fields 1^-6- 1 (ya) 
with the special value —b^^uji for their W3 momentum: 



/3n— 6 



^ ( n y-b-^^^iya)ll^aM) ) ■ (1.2) 
\ a=l 1=1 I 

The number of degenerate fields is of the order of 3n, which allows their worldsheet positions Ua 
to (approximately) correspond to the 3n components of the isospin variables xi - ■ ■ Xn- This will 
also allow fi„ to obey some differential equations which may be related to the KZ equations for 
Moreover, the tentative relation between il„ and (in will involve a simple twist function 

a<b i a i<j 

for some constants A, /i, z/ to be determined in terms of the level k of our 5^3 algebra; and the inte- 
gral transformation JC with integration kernel K{{xi}\{ya}, U\{zi}) which implements Sklyanin's 
separation of variables, and may therefore depend on the spins jj but not on the level k. We will 
then investigate the validity of the conjecture r2„ ~ /C • Qn^n = J dU dya K ■ Qn^n, or more 
explicitly 



an{{xi]\{z,}) I dUj{dya K{{x,]\{ya}M{z,]) ■ Qr^{{ya}\{z^mn{{ya]\{zi]) • (1.4) 



The meaning of the equivalence ^ here is that both sides obey the same differential equations. If 
true, this equivalence may then be promoted to a relation between physical correlation function of 
specific models, like the relation between the model and Liouville theory [|l|], but this is not 
the focus of the present article. This is why we do not worry about such details as the dependence 
of the correlation functions on antiholomorphic variables. 

The article will start with a brief review of the KZ equations and other Ward identities in 
conformal field theories with sIn symmetries, where we will explain how the Gaudin Hamiltonians 
appear in such equations. We will then review the KZ-BPZ relation in the 5^2 case; the reader is not 
advised to skip that section as the KZ-BPZ relation is presented in a form suitable for generalization 



to 5^3. In the 5^3 case, we will then find that the conjecture ( |l.4[ ) holds only in the critical level 
limit — >^ 3. 



2. Gaudin Hamiltonians in conformal field theory 

We will review how the Gaudin Hamiltonians appear in Ward identities obeyed by correlation 
functions in conformal field theories with an sIn symmetry algebra. The Ward identities associated 
to the stress-energy tensor T'^{z) lead to the KZ equations, which involve the ordinary Gaudin 
Hamiltonians. The Ward identities associated to the cubic field W'^{z) involve higher Gaudin 
Hamiltonians. 
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2.1 Knizhnik-Zamolodchikov equations 

The affine Lie algebra sIn is an infinite-dimensional extension of tlie simple Lie algebra s£j\f. The 
generators of si]\j, its structure constants Z"'', and its metric k"^ are defined by the relations 

[t'^^t^] = fff , k'^'' = Tr ft'' , fffl;'^ = 2Nk'"^ , (2.1) 

where here and in the following the trace Tr is taken in the fundamental representation, so that our 
metric k"'^ coincides with the renormalized Killing form of [^(13.13). The affine Lie algebra s£n 
can be formulated as the algebra of currents J^{z) with the operator product expansion 

r{z)j\w) = —. ^ + ff^^ + {rj'){w) + 0{z - w) , (2.2) 

[z — w)'^ z — w 

where the parameter k is called the level, and the normal-ordered product {J°-J^){w) is defined by 
the present formula. Conformal symmetry follows from the existence of a Virasoro algebra with 
central charge c = , generated by the Sugawara stress-energy tensor 

where J" J** is a shorthand for KabJ"'J^- The identification of T'^{z) with the generator of confor- 
mal transformations will be at the origin of the KZ equations. These equations are satisfied by any 
correlation function (1.1) of n affine primary fields ^^^{xi\zi) on the complex z-plane, where the 
spins ji label representations of s£n, the isospin variables Xi label the states in a given representa- 
tion, and the complex numbers Zi are positions on the Euclidean two-dimensional spacetime. The 
affine primary fields are defined by their operator product expansions with the currents J°'{z), 

r{z)^^ {x\w) = ' + Oil) , (2.4) 

z — w 

where D"^ provides a realization of the representation of spin j in terms of differential operators 
acting on the isospin variables x, so that [D"', D''] = f^'^D'^. We will keep this realization arbitrary, 
without committing to any particular choice of isospin variables. Let us however give an example 
of such a choice in the s^2 case: 

D- = — , D^ = x—-j , D+ = x^—-2jx. (2.5) 
ox ox ox 

The KZ equations are now obtained by inserting T"^ {z) into the correlation function il^, and us- 
ing the conformal Ward identity for T-^{z) on the one hand, and the affine Ward identities for 
( J" J")(z) on the other hand: 



T-^{z)Yl^^ 



Xi \ Zi j 
i=l 




2(k -N)^ z 

1=1 



(2.6) 
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where the subscript (z) in L)".^ indicates that it acts on the isospin variables Xj, and by definition 

P.(*) 



Lj[ is the p-th mode of T'^{z) acting on ^^^{xi\zi), according to 



Li<^>\x\z) = — Idwiw- zY^'^T\w)^Ux\z) . (2.7) 
^ 27rz .L 



Calling A'^ the eigenvalues of Lq, such that = A^.il„, we first deduce from eq. (2.6) the 



expression for Aj in terms of the quadratic Casimir C2{j) = D"'D"' of the sEn representation with 
spin j, 

A J _ C2{j) 



Now T'^{z) is assumed to generate conformal transformations, and in particular L;{^ (j)^n = 
(We define = as a derivative at fixed isospin variables.) Together with eq. 

(p^), this implies the KZ equations [ p^ 



{k-N)—nn = -H,nn , Hi = y^^^^, (2.9) 



OZi 



Zi 



The n commuting differential operators Hi are called the Gaudin Hamiltonians. Through its de- 
pendence on D^^^ and I?^^^ , each one of the n Hamiltonians involves all of the n isospin variables 
Xi, which makes the problem of their simultaneous diagonalization difficult. This difficulty will be 
solved by Sklyanin's separation of variables, which replaces the isospins Xj with new variables y^, 
and combines the Gaudin eigenvalue equations into an essentially equivalent set of equations, each 
of which involves only one of the new variables. 

2.2 Ward identities for the cubic field 

In addition to the quadratic invariant tensor k"^ = Tr t°-t^, it is possible to define the fully symmet- 
ric cubic invariant tensor 

^abc ^ rj^ ^^a^b^c ^ ^a^c^6^ _ ^^2.10) 

This tensor vanishes in the case of sl2, but not in the cases of sIn>z- It can then be used for 
constructing the invariant cubic field 

W\z)^\pdaUJyj''Jl){^) , (2.11) 

6 {k-N)2 



This generalizes the Sugawara construction (2.3), with however two substantial differences. First, 



while the field T-^{z) is interpreted as the generator of conformal transformations, there is no such 
geometrical interpretation for W'-'{z). Second, while the field T-^{z) obeys a Virasoro algebra. 



the field W'^{z) does not obey the higher W3 algebra [11]. In other words, while the Virasoro 
algebra can be realized as either a coset of 5^2 or a subalgebra of the enveloping algebra of s£n>2 
(albeit with differing central charges), the VVs algebra is a coset of 5^3 but not a subalgebra of the 
enveloping algebra of s^Ar>3. 
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In analogy with eq. (2.6) we now have 



i=l 




+ 



+ 



2,W 



{z - Zj)3 {z - ZiY Z 



Zi 



Z — Zi V — Vn ^-^ - 



1 = 1 



= 1 



m=l 



Z 



where by definition is the p-th mode of W'\z') acting on $-'»(xj[zj), according to 



W^^\x\z) 



1 
27ri 



dw {w - zY"^'^W\w)^\x\z) 



(2.12) 



(2.13) 



Calling q'^ the eigenvalues of Wq , such that ^^j^p^n = Qj^^n, we first deduce from eq. ( |2.12 ) the 
expression for q-J in terms of the cubic Casimir Csij) = dabdD" D'' D" + D^'D^D^) of the s£n 
representation with spin j, 

\pC3U). (2.14) 



6'^ 



We further deduce 



1 
2 

\pH'lnn 



(2.15) 
(2.16) 



where the differential operators H[ and ff -' are higher Gaudin Hamiltonians, whose explicit expres- 
sions in terms of Z)^^^ can easily be derived from eq. ( 2.12 ). But, in contrast to L^_i, the operators 
W:!_i and W:!_2 iiot interpreted as differential operators with respect to z. The equations ( |2.15 ) 
and ( 2.16 ), which generalize the KZ equations, are therefore not differential equations, and they 
will therefore not help us test our conjecture. Nevertheless, they will naturally appear in certain 
formulas. 



3. Review of the 5^2 case 

In this section we will review the relation between the s^2 KZ equations and BPZ equations. This 
was originally found by Feigin, Frenkel and Stoyanovsky [||], using Sklyanin's separation of vari- 
ables for the s^2 Gaudin model [^. However, the original derivation relied on a particular choice 
of the isospin variables. This choice of isospin variables makes the result remarkably simple, but 
has no analog in the §£3 case, as we will show. We will therefore reanalyze the 5^2 case, using 
whenever possible objects which do have analogs in the 5^3 or even sIn cases. We will present 
systematic derivations of their relevant properties, which will help clarify whether and how they 
can be generalized to the 5^3 case. 

3.1 Separation of variables for the s(.2 Gaudin model 

Let us consider a system of n representations of 5^2 with spins ji ■ ■ ■ jn- Consider the associated 
quantum variables D'^.^ such that = SijffD"^^ with D'l^^D''^^^ = C2{ji). The system 

comes with parameters zi • • • Zn- Sklyanin's separation of variables for this system involves three 
ingredients: 
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1. A function B{u) of an arbitrary variable u (the spectral parameter), whose zeroes are the 
separated variables yi, so that B{yi) = 0; 

2. Another function A{u) such that pi = A{yi) is the conjugate momenta to y,; 

3. Akinematical identity, called the characteristic equation, which for any given i relates yi and 
Pi- 

We now briefly review the construction of these three objects in the s^2 case. They are built from 
the s£2 Lax matrix 



(i) 



i=l 



U - Zi 



(3.1) 



whose matrix elements Jq (n) obey the identity 

(n - v)\ll[u)J%{v)\ = 5ir<{u) - 5lUu) - S^I^iv) + 6}r^iv) . (3.2) 
With the particular choice eq. ( |2.5| ) for the si2 isospin variable x, the si2 Lax matrix is explicitly 
' iV" -J— (^.^ - -i-'] V'' -J^ ^ 



I(u) 



(3.3) 



Now choosing 



it is easy to check that 



[B{u),B{v)]=0 , [A{u),A{v)]=0 
{u - v)[A{u), B{v)] = B{v) - B{u) . 



(3.4) 

(3.5) 
(3.6) 



These relations ensure that the operators yi defined as the zeroes of B{u), and pi = A{yi), do 
satisfy 



[yj,yj] = o , \pi,yj\ = 5ij , [Pi,Pj] = o. (3.7) 

In particular, \pi,B{v)] = agrees with B{v) oc [|'()^_^'.) ■ There is however a problem 

of operator ordering in the expressions A{yi) and B{yi), because the separated variables yi are 
operators. This problem is dealt with in reference [^. We will ignore it in the forthcoming heuristic 
derivation of the characteristic equation. Let us start with det {A{yi)\d — I{yi)) = 0, where id is 
the identity matrix. (The determinant of a matrix whose first line vanishes is zero.) This implies 

P' 



-^^ — ^{Ialp){yi) = 0. This characteristic equation can easily be rewritten as 




(3.8) 



where is of course a Gaudin Hamiltonian ([2.9|), and C2{j) is the quadratic Casimir of a spin-j 
representation. 
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Functional space interpretation. We now wish to consider the quantum variables D^^^ as dif- 
ferential operators acting on functions ^({xj}) of isospin variables Xj. (An example of such a 
realization was given in eq. (|2.5[).) Similarly, the separated variables and their associated mo- 



menta pi may act on functions ^{{yi}), in particular = The separation of variables 

{xi} — > {yi}, U (where the extra variable U will be defined shortly) is then intepreted as an integral 
transformation K, such that 

^{{xi})=}C^{{ye},U) = jdU I l{dye K{{xi}\{ye},U) ^{{ye},U) , (3.9) 

where the kernel K is characterized as a common eigenvector of the commuting operators B{u) 

B{u) - C/MiZL^'j K{{x,]\{y,], U) = . (3.10) 
The simultaneous diagonalization of the Gaudin Hamiltonians Hj, namely the set of equations 



{Hi — Ei)^ = 0, can now be reformulated using the characteristic equation (3.8), which implies 



f 1 V- C2(je) ^ Ee \ ~ 

hr^-oEr^^-E — ^^ = 0, (3.11) 

The solutions of this equation can be found in factorized form ^' = Hi V'(yi)- This justifies the 
name "separation of variables" attributed to the change of variables Xi —>■ yi. 

Some remarks. Finding the kernel K by the simultaneous diagonalization of the operators B{u) 
is easy in the sl2 case because B{u) = (n) is a sum of n commuting operators, so that we have 



K{{xi}\{y(}, U) = nr=i ^) where the the equation on ki is obtained from eq. ( |3.10D 

in the Umit u — > zi. 

{nlDl^+^l^h{x,\{yl},U) = Q , EE [/iM^i^^ . (3.12) 



For example, if the isospin variables are chosen as in eq. (2.5), then we find ki = e This 



suggests that we could use other isospin variables /xj such that D'^^^{t'^)\ = —fn, then we would 



find ki oc 6{fli — fii), so that we could explicitly perform the integrals in eq. ( |3.9| ). This would 
lead to ^'({/ij}) oc ^({y^},?/) with simple proportionality factors, as the change of variables 
{fj-i} —>■ {{ye}, U) would now be local and described by the functions /Uj({y^}, U). More generally, 
for any choice of isospin variables, the kernel K will be of the type 

n 

K{{xi]\{yi], U\{zj}) = J{h {xi\ {/X,}) , (3.13) 



where /Xjdy^}, C/|{zj}) is defined in eq. (3.12), and we made the -dependence explicit. Thus, in 
the case, the kernel K can be determined explicitly, and this is because the operator B{u) is a 
Unear function of the Lax matrix I{u). 

Let us finally be more precise about the number of variables y^. They are defined as the zeroes 
of a rational function B{u) which, barring extra constraints, has n poles and degree —1. Therefore 
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we must have n — I such variables, and the nth variable U is the eigenvalue of — Y17=i ^\ 



In conformal field theory applications, we however impose the extra constraint X^iLi 



a 

0, so 



that B{u) has degree —2. This yields n — 2 variables {ye}e=i---n-2, and U is the eigenvalue of 
3.2 The s£2 Knizhnik-Zamolodchikov equations in Sklyanin variables 

We just saw that Sklyanin's separation of variables is useful tool for simultaneously diagonalizing 
the s£2 Gaudin Hamiltonians. This problem is closely related to the problem of solving the KZ 
equations ( [2.9[ ), which are obtained by replacing the eigenvalues of the Gaudin Hamiltonians Hi 
with — {k — 2)j-. This suggests that it may be interesting to rewrite the KZ equations in terms of 



Sklyanin's variables. To do this, we will use the characteritic equation (3.8) which such variables 



obey, and apply it to /C^^$7„, which is a function of {yi}, so that pilC~^Qn = ^^~^^n- While 
itself just a kinematical identity, the characteristic equation then allows us to reorganize the KZ 
equations as 




(3.14) 



where we drop the index from y^, and we use Aj = — 2{k-2) fro^n eq. ( [2.8[ ). We still have to 
perform the change of variables on the ^^-derivatives at fixed isospins, i.e. to rewrite /C~^^/C in 

. This is rather easy because of the particular form of the kernel ( |3.13| ), where 

the dependences on {ya}-, U and {zg} are channeled through the particular functions {/ij}. This 
implies that the integral transformation ( ^!9| ) just adds first-order differential operators to 



terms of 

dz, 



a 

dzi 



8zt 



, so that 



6z£ dz£ ^ dz£ 



d dU 

dya dzi 



d_ 

dU 



(3.15) 



Denoting {ya} = {y, {yb}}, we obtain the KZ equations in Sklyanin variables, 



k-2dy^ 



^ 1 / d d_\ ^ 1 r d d 

^ y — Zi \dzi dy ' ^ k — \ ^H'^ 



y-yb \dyb dy 



Jt 



K,-^ a„ = . (3.16) 



r{{y- zif / 

In this equation the variables are no longer separated, as the variables ?/b appear in addition to y. 



3.3 Comparison with Virasoro null-vector equations 

In the previous subsection, we have studied the KZ equations in a CFT with an sl2 symmetry 
algebra at level k. We will now compare them with null-vector equations in a CFT with a Virasoro 
symmetry algebra at central charge c = 1 + 6(6 + 6^^)^ where 6^ = This is the Virasoro 
algebra which would be obtained from our s^2 algebra by quantum Hamiltonian reduction (see for 



-9- 



instance [|n|]), although that reduction does not explain the relation between differential equations 
which we are about to review. 

The Virasoro algebra can be formulated in terms of the stress-energy tensor T{z), which obeys 

1 „ 



Tiz)T{w) 



+ 



2T(w) 



+ ^ + 01 



w 



w 



w 



(3.17) 



Primary fields Va{w] 
by 



of momentum a and conformal dimention 



a 



w 



W 



> + b ^ — a) are defined 



(3.18) 



{z — wY 

This definition does not distinguish the primary fields Va and 1^6+6-i_q,, which have the same 
conformal dimension. These fields are therefore assumed to be proportional, with a proportionality 
constant called the reflection coefficient. This Z2 symmetry can be understood as the action of the 
Weyl group of SI2 on the space of the momenta a. 

The Virasoro representation generated by the degenerate field V_i_ is known to have a nuU- 
vector at level two Namely, (L_2 + = 0, where the modes Lp are defined as in eq. 



(2.7). This implies that correlation functions involving such a degenerate field obey the Belavin- 



Polyakov-Zamolodchikov equation 012| ] 

- 1 Q 



+E 



y 



n d 



+ 



n 



A, 



y.^{y)^yaM)) = (3.19) 



Curiously, this equation is formally identical to the variable-separated KZ equation (3.14). The 



meaning of this formal similarity is not clear to us. The KZ equations in Sklyanin variables ( |3.16| ) 
actually involve n 
tions of the type 



2 variables yi • • • yn-2, therefore we should rather consider correlation func- 




= ( n ^-A.(ya)li^a.(^») ) . (3.20) 

\a=l i=l 

We then expect such correlation functions to be related to il„ ( [Ell ) as in equation (L4). That 
equation means that the twisted BPZ equations satisfied by Qn^n are identical to the KZ equations 



in Sklyanin variables ( [3.16| ). This can indeed be checked by explicit calculation, provided we 
correctly specify the function ©„ as well as the relation between 5^2 spins ji and Virasoro momenta 
CKj. Requiring that the a—j relation is compatible with the respective Weyl symmetries j —>■ —j — l 
and a b + b"^ — a, and that conformal dimensions Aj - 
by a constant shift, determines the relation 

J_ A - A-^ i ^ 

26 ' - + 2 + 462 



''^^^2^ eq. ( |2.8[ ) and Aq, are related 



a 



Kj + 1) + 



(3.21) 



We still have to specify the values of the parameters A, /i, u in the ansatz (1.3) for the function 0„. 



We could determine these values by requiring the twisted BPZ equations to agree with eq. (3.16), 
and we would find 

1 _ 1 _ 1 

262 ' ^--^ ' 



A 



(3.22) 
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There are simple concurring arguments for the values of A and u. First, the value of A is determined 
by the requirement of continuity of QrS^n at Ua = Uh- This requirement plays an important role in 



the boundary model [^]. Second, the value of u follows from checking equation ( |1.4| ) in the 
simplest case n = 2, when there are no ya variables and no BPZ equations. 

Let us now comment on this twist function 0„ and its relation to free field correlation func- 
tions. In this paragraph we will consider full correlation functions with dependences on both holo- 
morphic and antiholomorphic variables, and the full twist factor which is thus |0nP- With the 



above values ( 3.22 ) for X,fi,u, we observe that the inverse twist factor coincides with 

the free field correlation function formally obtained from 0,^ by taking the fields Vai{zi) to have 
momenta Q-i = ^ instead of = 6(jj + 1) + ^. This means 

/n-2 n \ 

\®n\-^ = l'Y{V_i_{ya)J{VAzi)) ■ (3.23) 

\a=l i=l I 



This interpretation of 0„ plays a role in a recent proof of the FZZ conjecture ||13|], see also [|1^ 
For now, let us explain the origin of this observed relation by studying the 6 — > limit of the 
iJ^-Liouville relation. This relation can be written as Q.n. ~ lCiC\Qn\'^^n > whose factors we now 
analyze: 

• The Liouville correlation function f2„ reduces to (na=i ^-^iVa) IlILi ^J-(^i) / as 6 — > 
0. And it turns out that this coincides with a free field correlation function, because the 
momentum conservation condition is obeyed. Namely, the sum of the momenta is (n — 2) x 
~ \ which coincides with the dominant term in the Liouville background 
charge \ ^h. Therefore, according to standard path-integral reasoning in Liouville theory 

free ^ 



Q, we have 17„ ~ ( Da-i ^-^ (^a) 11^1 ^-^ (^i) > =^ 



n(;/a-z») 



n(ya-j/6)n(zi-^3) 



where i?„ is 6-independent. 



The i^g correlation function r2„ is expected to have a finite "minisuperspace" limit [16] as 
6^0 which is equivalent to ^ oo where k is the level. 

The separation of variables K is 6-independent by definition. 

So the twist factor |0nP must absorb the 6^0 divergence of the Liouville correlation func- 



tion which implies the relation ( |3.23| ) and the values ( [3.22| ) for the parameters A, ^, v. 
(This reasoning does not exclude the presence of extra terms in A, /i, v which would be finite 
in the 6 — > limit.) 

This concludes our reminder of the KZ-BPZ relation in the sl2 case. In the next section we 
will analyze the siz KZ equations along the same lines. 



4. The case 

4.1 Separation of variables for the s^s Gaudin model 

To the best of our knowledge, the full quantum separation of variables for the 5^3 Gaudin model has 
not been derived yet. By the full separation of variables we mean the determination of A(u), B{u) 
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and a characteristic equation, like in the 5^2 case. ' Sklyanin did however derive the full separation 



of variables for the classical 5^3 Gaudin model [18]. In order to derive the quantum version, we 



will use Sklyanin's separation of variables for models with an si^ Yangian symmetry m\, see also 



[19] for a generalization to s£]\[. This Yangian symmetry is present in the Gaudin model, which 



will allow us to derive its quantum characteristic equation from the Yangian's. 

Yangian symmetry. As in the 5^2 case, the variables of the s£n Gaudin model can be com- 



bined into an s^at Lax matrix I{u) (p.lj) obeying the relation (3.2). It is however possible to 



combine the variables into another sIn matrix, which depends on an extra parameter r], 

= id + r,I{u) + ^r?2 : /2 : (u) + ^r]^ : : (u) + ■ ■ ■ , (4.2) 
2 b 

where the definition of the normal ordering in : : (u) and : : (n) follows from the chosen 
ordering of the factors of Y{u). This object can be shown to obey the Yangian algebra 

(n - v)Y^{u)Y^iv) + vY^{u)Y^{v) = {u - v)Y^{v)Y^iu) + vY^{v)Yjiu) . (4.3) 

Sklyanin's separated variables m for the Yangian [||] are defined as the zeroes of a function 

B^{u) = Yi{u)Yi{u)Yi{u - 7?) - Yi{u)Yi{u)Yi{u - r?) 

+ Yi{u)Yi{u)Yl{u - 7?) - Yi{u)Y^{u)Yi{u - r?) , (4.4) 

while the conjugate variables are given by Xi = [yi) where 

A^{u) = Yliu) - Yiiu - vr%\u - rj)Y,\u) , (4.5) 

Let us point out that interesting structural insight into these formulas for {u) and B^{u) was 

obtained in [pO|], based on general properties of matrices with non-commuting elements. The 
functions A^^ (n) and (n) obey the commutation relations 

[A''iu),A^iv)]=0 , [B^{u),B^{v)]=0 , ^^[A^ (u), B^" (v)] 

V 

= B^{u)A^{v) Yiiu - 7^r'Yi{u)-%\v - r^)Yi{v) - B^iv)A^{u) , (4.6) 

so that 

[y^,yj]=0 , [Xi,y,] = -r,6i,Xi , [Xi,Xj]=0. (4.7) 
The quantum characteristic equation is then 

Xf - Xfh{y,) + XMy^ -r])- d{yi - 2r?) = , (4.8) 



'a different approach was proposed in ||l7|], which consists in trying to use the SI2 separation of variables in the sl-j, 
case. This approach requires a particular choice of isospin variables. The results are complicated. 
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with the invariant operators ti(n), t2{u) and d{u) defined as ^ 

h{u)=TTY{u) , t2{u)=TvY{u) , d{u)6Z=Yl^{u)Yj{u + r]) , (4.9) 

where the matrix Y is constructed by transposing the quantum comatrix of Y. For instance, 
I3 (u) = —Y^{u)Y^{u + r]) + Y^{u)Y^{u + tj), where the /^-shifts are the manifestation of 
the quantum character of the comatrix whose 2 matrix element we just wrote. Operator ordering 
issues in expressions like t2 {ui — rj) are resolved by inserting the operator yi from the left. 

From the Yangian to the Gaudin model. We will now construct objects A{u), B{u) and a 

quantum characteristic equation for the 5^3 Gaudin model. Such r^-independent functions of the 
matrix I{u) will be obtained by expanding the con^esponding objects for the si^ Yangian algebra 
in powers of rj. We find 

^^(n) = l-r?A(n) + 0(7/2) , ^(u) = _/i + l^L , (4.10) 

-'3 

B^{u)= r^^B{u) + Oiv^), B{u)=llllll-Illlll + Illll\-Illlll, (4.11) 

where we omitted the spectral parameter u in Ia{u), and we point out that our formula for A{u) is 
free of ordering ambiguities because /|(it) commutes with both lf{u) and /|(n) . The commutation 



relations (4.6) for A (u) and B (u) imply the analogous relations 



[A{u),A{v)]=0 , [Biu),Biv)]=0, (4.12) 
{n - v)[Aiu),Biv)] = B{v) - , (4.13) 

-'3 V^J^3 



which may be compared to the corresponding relations in the SI2 case eq. (3.6) 



Let us rewrite the characteristic equation (18) as: 



{Xi - if - {Xi - if [ti{yi) - 3] + {X, - 1) [t2{y^ - r]) - 2ti{yi) + 3] 

+ [1 - hiyi) + t2{y^ - 7?) - d{yi - 2??)] = . (4.14) 

The leading behaviour of this equation as 77 ^ will turn out to be 0{rj^). To compute this 
behaviour, we of course need to compute the behaviours of Xi and yj as 77 — > 0. It turns out that we 
only need the 0{r]) behaviour of X^. We therefore define the variable pj by = 1 — r]pi + 0{rf). 
As for yi we only need need the leading 0(1) behaviour. To this leading order, the zeroes of B^{u) 
coincide with those of B{u), so that we do not need distinct notations and call them all yi. The 
most complicated part of the calculation however does not involve such subtleties, but rather deals 
with the last term in eq. ( |4.14D , 

1 - ti(n) + t2{u -T])- d{u - 2ri) 

= (1 - Yliu - 2r,)){Yi{u -v)- l)Yi{u) + {Yl{u - r,) - Yl{u - 2j^))Yi{u) 
+(y/(u - 2r/) - l)Yi{u - r])Yiiu) + - 7?) - l)Yiiu) + (1 - YHu)) 

+iYi{u - 2r?) - l)Y,^iu - v)Y^{u) - Yi{u - 2^)Yf{u - r^)Y^{u) 
-Y^{u - 2r])Y^{u - ri)Yi{u) + Yf{u - 2r])Yi{u - v)yiiu) - Y^iu - v)Yiiu) 

= [-iiiiii + iiim + ipM + mil - iiiiii - iim . 



-13- 



where we omitted the spectral parameter u in Ia{u), and used the s^s-defining relation l| + /| + 
Ig = 0. We then obtain the following quantum characteristic equation of the 5^3 Gaudin model: 



v^-P^■\{I^Ip){yi) + -AI^Ip)'{y^ 



6 



(4.16) 



Notice that the particular cubic invariant which appears in this formula is related to the fully sym- 



metric invariant tensor dabc eq. (|2.10| ). Using the definition ( p.lt ) of I{u), we indeed have 



I^Ip^+ip^i:){u) = -d. 



^abc 



1=1 



7/ _ r.n ^-^ 



(m) 



£=l m=l 



(4.17) 



This could further be expressed in terms of the higher Gaudin Hamiltonians of Section [2.2[ , so that 
the characteristic equation could help simultaneously diagonalize these Hamiltonians. 

Some remarks. Like in the s£2 case, Sklyanin's change of variables can be interpreted as an 



integral transformation /C (3.9) acting on a functional space. The kernel K of K, now obeys 



B{u) - U 



u ■ 



ye) 



Uiiu-z^y' 



K{{x^}\{yi},U)=0 



(4.18) 



However, the simultaneous diagonalization of the commuting operators B{u) is now a difficult 
problem, as B{u) is now cubic and not linear in I{u), and thus no longer a sum of n commuting 



operators. Therefore, the kernel K is no longer of the form ( 3.13 ). Certainly, no choice of isospin 
variables exists such that the kernel K has a simple expression. Another difference with the s£2 case 
is the counting of variables: generic functions of the si^ isospin coordinates Xj should correspond to 
functions of not only yi and U, but also of two extra variables. These extra variables are necessary 
for the transformation /C to be invertible. We will neglect this issue^, as well as the issue of precisely 
defining the relevant functional spaces, and we will assume /C to be invertible. 

Let us finally determine the number of separated variables - that is, the number of zeroes 
of B{u). Barring extra constraints, this is of course 3n — 3. In conformal field theory applications, 
we however impose the extra constraints Y17=i ^(i) ~ 0' ^^^^ ^(^) degree —2. This does 
not immediately imply that B{u) (eq. ( |4-.11| )), which is cubic in I{u), has degree —6, because 
Y17=i ^{i) ~ only holds when directly applied to a physical correlation function, and the matrix 
elements of I{u) generically do not commute with each other. Rather, the degree of B{u) depends 
on its precise form and should be evaluated by explicit calculation. We find that each one of the four 
terms of B{u) has degree —5, while B{u) itself has degree —6. This means that there are 3n — 6 
separated variables. Therefore, as in the s^2 case, the number of separated variables vanishes for 
n = 2. 

4.2 The sis Knizhnik-Zamolodchikov equations in Sklyanin variables 

Let us consider a conformal field theory with an si^ symmetry algebra. The Ward identities consist 



in the n KZ differential equations (2.9), plus 2n extra non-differential relations (2.15) and (2.16), 



which express ^--^ and W_2 in terms of differential operators acting on isospin variables. Let 



-2,» 



"A construction of the extra variables seems to be available in the article ili 
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us reorganize all these relations by injecting them into the characteristic equation of the quantum 



5^3 Gaudin model (4.16). The result is schematically of the form 



+ 



3)^ • T-^(y) - }:{k - 3)dT\y) - -W\y) 
dy 2 p ^ 



(4.19) 



where the constant p was defined in eq. (2. 1 1 ). Explicitly, 



dy'^ 



+ {k-2,) 



d_ 

dy 



i=i \ 
+ -ik 



— /c-i— /C + — 

Zi OZi [y 



3) 



'/C 



1=1 



{y 



- ZiY^ 5zi^^ {y 



2Ai 



y 



Zi 



+ 



{y 



+ 



4 



{y - Zif 



JC-^Qn = , (4.20) 



where Qn is still an n-point function of the type 

In this equation, the terms involving ^ ^.-^ and W;^^ (j) refer to correlation functions involv- 
ing descendents of the primary fields ^^p\z). We have little control over such non-differential 
terms, and we would like to ignore them in the following. This could be done by considering 
appropriate linear combinations of our 3n — 6 equations. (Remember that the variable y spans 
the 3n — 6 separated variables {ya})- We will for simplicity adopt the alternative approach of 
working modulo the unwanted terms. Let us make this precise by defining the space Vs of differ- 
ential operators in ya, zi (including functions of ya, Zi) which are symmetric under permutations of 
{yij 2/2 • • • ysn-e}- For any choice {ya} = {y, yt} of a distinguished variable y we further define 



^2(!/) 



1 



y 



-Vs + 



1 



1 (y - Zi)' 



(4.21) 



By a simple counting of variables it can be realized that any differential operator which is sym- 
metric under permuations of {yb} does belong to J^3{y) = Ya=i ip^'^s + E"=i (.y\)'i '^s + 



En 
1=1 



Vg. But it does not always belong to J^2(y)> so we can define a nontrivial equivalence 



as the equaUty modulo T2{y). Thus, equation ( 4.20] ) simplifies to 




(4.22) 

Having thus eliminated W;[^ r--, and 1^^2 (i)' ^^^^ ^^^^ operators which we recall are 



-derivatives at fixed isospin variables. We expect K. ^-r-K, to be a combination of the operators 



and although we do not know how to compute it. And it is not clear whether K. ^j^K, 
is a first-order differential operator, as happened in the case (see eq. (p.l5[)). Nevertheless, we 



do know that /C ^^f-AI is independent from the level k, which is a parameter of our conformal 
field theory but neither of the Gaudin model nor of its separation of variables. Therefore, we will 



still be able to extract useful information from eq. ( |4-.22[ ), a sum of terms with various power-like 
dependences on {k — 3), by considering all terms which are not linear in {k — 3). 
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4.3 Ws null-vector equations 

Let us first briefly explain why we try to relate conformal field theories with an 5^3 symmetry at 
level k to theories with a Wa symmetry at central charge c = 2 + 24(6 + 6^^)^ where 



1 



/c-3 



(4.23) 



A theory with an sl^ symmetry like the 5^3 (R) WZW model can be written in terms of eight 
quantum fields, as 5^3 is eight-dimensional. However, affine sl^ highest-weight representations 
are parametrized by just two numbers, namely the two components of the 5^3 spin j. This suggests 
that the non-trivial dynamics of the theory really take place in a two-dimensional space, where j 
would play the role of the momentum. There exists such an s^3-based theory which involves just 
two interacting quantum fields: the conformal 5^3 Toda theory, which has a symmetry algebra. 
The correct parameter h for this algebra is suggested by the Drinfeld-Sokolov reduction, which 
realizes yV3 as a kind of coset of the sl^ algebra. 

algebra and primary fields. Referring to the review article [ jlTl ] for more details, we recall 
that the W3 algebra is spanned by the modes of the fields T{z) = J2nez -^n-z ""^ and W{z) = 
Ylin&'L WnZ~'^~'^. Let US Write the defining relations of the W3 algebra in the form of commutation 
relations for the modes L„, Wn rather than operator product expansions for the fields T{z),W{z), 
as this form is more convenient for finding null vectors in representations: 

[Lm, Ln] = (m - n)Lm+n + -^m{w? - l)5m+nfl , (4.24) 

[L^, Wn] = (2m - n)Wn,+n , (4.25) 

(22 + 5c) m — n 2 2 1 a\ t 1/ na ,a 

H 15~*-"^ + n - ^mn - 4)Lm+„ + -{m - n)Am+n , (4.26) 

where we introduce, using the normal ordering : LmLn '■= LmLn ifm <n, 

Am = y^-- LnLm-n ■ +^XmLm with \ ^'^^ ^1 r t ' (^.27) 

^ 5 [x2^+i = (^ + 2)(1-^) 

A primary fields Va of the W3 algebra of momentum a, conformal dimension and charge qa is 



defined by its operator product expansions with T{z) eq. (3.18) and W{z): 



W{z)Va(w) = - ^ + -7 yT ^ ^C'(l). (4.28) 

[z — wy [Z — wY Z — W 

The momenta a now belong to the two-dimensional root space of the Lie algebra 5^3. A basis of 
this space is provided by the simple roots ei, 62 whose scalar products appear in the Cartan matrix 

( {e2'el) (62 62) ) ~ ( -1 )• '^^y ^1^° basis wi = |ei + |e2, u)2 = \ei + |e2 

such that {ei,u}j) = 6ij. We decompose the momenta along this dual basis: a = aitoi + a2U!2, 



-16- 



and we introduce the vector Q = {b + b ^)(ei + 62). The conformal dimension and charge are 
parametrized in terms of the momentum as 



^{a,2Q-a) 



27 



[ai -a2][2ai + 02 -3(6 + + 2^2 -3(6 + 6-1)] ^ 



(4.29) 
(4.30) 



W3 degenerate fields. Let us now justify the choice of the field V_ 



in the correlator $7^ 



(1.2) which appears in our conjecture. We wish Qn to obey third-order differential equations, 
which would correspond to the sis KZ equations in Sklyanin variables. This suggests that we use 
the simplest non-trivial degenerate fields, which have null vectors at levels 1, 2 and 3. But there 
are actually four such degenerate fields, with a € {—bcoi, —buj2, —b~^uji, —b''^uj2}, whereas we 
want only one of them to appear in 0.^, because the original isospin variables are invariant under 
permutations of the Sklyanin variables. 

By analogy with the 5^2 case, we focus on the fields and V_i,-i^^^, whose momenta go 

to zero in the critical level limit /c ^ 3. They are related to the other two fields by the Ws algebra 
self-duality b b^^, which is however not an invariance of the 5^3 algebra. And they are related to 
each other by the Dynkin diagram automorphism uji <-> uj2 of 5^3, which acts on general primary 
fields Va as (Aq,, qa) (Aq,, —Qa)- This symmetry does have a counterpart in the separation of 
variables for the 5^3 Gaudin model. The construction of the separated variables was indeed based 
on the introduction of an si^ Lax matrix I{u) ( |3.l| ), so that 5^3 generators act in the fundamental 
representation. But we could alternatively have used the antifundamental representation, which is 
related to the fundamental by the Dynkin diagram automorphism. With our conventions, our choice 
of the fundamental representation will turn out to correspond to the choice of the degenerate field 
V_i,-i^ of the W3 algebra. The three corresponding null- vector equations are [p2|] 



[^^-1 + (I + ^-1] ^- 
[iW.2-iL.2-bL ' 

[iW-3-{'^ + i)L. 



6-1^1 



IJ ^-b--iJi 

3 + 6L_iL_2 



0, 

= , 



b'Ll,] 







(4.31) 
(4.32) 
(4.33) 



The last null-vector equation implies that any correlation function with one degenerate field obeys 

El (Vl,-i^,(y) nr=i y»M)) = 0, where 



262 ^ 664 / ^ 



i=l 



1 



d 



+ 



2A, 



(y - ZiY dzi {y - Zi)^ 



+ 



63 



E 

i=l 



+ 



+ 



Zi {y - ZiY [y - Zif 



(4.34) 



This may be compared with eq. ( |4.20| ), which is formally similar, or even identical if the term with 
coefficient is absorbed into the other terms by redefining VF„i (j) and q^^. Like in the 5^2 case, 
the meaning of this formal similarity is not clear. 
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Now the equations obeyed by correlation functions with several degenerate fields like 17„ 



eq. (1.2) are significantly more complicated than Ei, because eliminating W-i, W-2 descendents 



of the degenerate fields requires the use of the first two null-vector equations ( |4.31|J4.32D . Still 
denoting {ya} = {y, yb}, we obtain the equation -E2^n = with 



E2 = El + 



-y 



y 



1 d 

+ 



Vbdyf ' b'^dy ^ 



E 



d 



+ 



y-Vb dyb {y - Vb) 



2 1 

0,1 



364 



E 



1 



(y - Vb) 



d_ 

dzi 

_d_ d_ 

dyb dy 



+ 



yb 



+ 



+ — V 1 / d ^ 

f^^(y- yb)(yb - yc) \dyc 



1 1 

62 ^64 



A 



-b- 



+ -^^q-b- 



^ {y 

b 



yb-y 
1 



yb) 



(4.35) 



where 



-1 



4 

362 



2763 



(4 + 362)(5 + 362 



Relating W3 momenta to 5^3 spins. 



(4.36) 
with the KZ 



In order to compare the equation E2^n 
equations in Sklyanin variables ( 4.22 ), we should specify how we relate si^ primary fields <I>-' (/i|z) 
to W3 primary fields Va{z). We are looking for a relation between a and j which translates into a 
simple relation between (A^, qa )and (Af,g/). We propose 




(4.37) 



where we use the following expressions for (A^, qj ) defined in eqs. (18) and ( 2.14 ) 



A:^ 



11. 

;-(j,j + 2ei + 2e2) 



A;-32 

-^-^T^bl - j2] [2(jl + 1) + (i2 + 1)] [(jl + 1) + 2(i2 + 1)] , 



{k-sy. 



27' 



(4.38) 
(4.39) 



where the components (ji , j^) of the spin j are defined as j = jiuJi + j2'^2- Notice that our relation 
between a and j maps the principal unitary series of 5^3 representations j G — ei — 62 + to 
the representations which appear in the physical spectrum of conformal si^ Toda theory 
a S Q + iM2. Such choices of a or j lead to real values of (A, g) if A; > 3. 

However, there does not need to be any relation between the 5^3 creation operators W^^ , W;l2 
and their Ws counterparts W^i,W-2- While relating L'/^i = ^ to L_i = though difficult 
in practice, is in principle a simple matter of performing the change of variables, there is appar- 
ently no principle which would determine how W^i , would behave through the change of 
variables. This is why we work modulo J^2{y), ignoring the non-differential terms which involve 
such operators, and being left with differential equations. Now the presence of degenerate fields in 
correlation functions of W3 fields does not necessarily lead to differential equations, a fact which 
makes conformal 5^3 Toda theory much more complicated than Liouville theory [0]. Differential 
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equations actually appear provided the number of degenerate fields is large enough. We are insert- 
ing 3n — 6 degenerate fields together with the n generic fields Fq-, which is enough for 
eliminating the 2n terms W_2,{i) and being left with n — 6 differential equations. 

Twisting Ws null-vector equations. Finally, we should determine the twist factor 0„ which 



appears in the conjecture ( |1.4| ), so as to be able to compute 

E3 = enE2Qn^ such that ^3 • e„fi„ = 



(4.40) 



The values of the parameters A, v can be derived as in the si2 case. Requiring continuity of @n^r, 



at Va = Ub implies A = 2A_^-ij^^ — A_2b-iaji = and requiring that the conjecture (L4) holds 



in the case n = 2 implies u = 2Aq — 2A^ = ^ + 4, see eq. ( [4.37 ). Notice however that this 
only determines v up to 6-independent terms, as the unknown 6-independent kernel K may also 
contribute. 

These constraints leave the parameter jj, arbitrary. We will obtain an ansatz for fi, and confirm 



the values of A and u, by generalizing the relation ( [3.23[ ) between Gn and free field correlation 
functions which was observed in the s£2 case. In the 5^3 case the analogous relation is 



/3n— 6 



free 



This ansatz leads to the values 



n y^b-^^Aya)ll^b 



b-'>-{ei+e2){^i) 



(4.41) 



a=l 



1=1 



2 

362 



2 
62 



(4.42) 



These values will turn out to be the only ones such that, modulo T2{y), the only non-differential 



terms in are of the type 



i . This is a rather non-trivial requirement as many non-differential 



terms can potentially appear (cf Appendix |A.1| ). Working modulo ^2(2/) eq- ( |4-21| ), and using the 
relation ( |4.37 ) between 5^3 and W3 representation data, we indeed compute 



-5/3 



1 ^ 1 1 ^-^ 

9^+52^2 + ^1^1 + ^2^ 



n n i J 1 aJ 



(4.43) 



where we introduced two differential operators Di and D2 of respective orders 1 and 2, which 
depend neither on the field momenta aj nor on the model parameter 6, 

/ \ 2 

d 



V ^ 2 1 V ^ 1 

^-^ (y — 2:j)2 dy > ^-^ 11 — ' f^ii 



f y- Zi dy 



-^E Yl 7^ - TT - 7^ - TT ^4.44) 

^-^ y — Zi ^—^ y ~ Ub \dyb dy ' ^-^ n — — n- \ '^■"i 'in 



d 



1 d f d 
d 



y -ybVb-yc \dyb dy 



bj^c 

d d 



^ y - Zidyi^dzi^^dyJ ' ^ y - yb \dyb dy J \dyb ' ~ dy 



d d 
+ 2- 



(y - ybY dy 



(4.45) 
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4.4 Comparing 5^3 Knizhnik-Zamolodchikov equations with W3 null-vector equations 



We are now in a position to test the conjecture (1.4) by comparing the KZ equations in Sklyanin 



variables (4.22), which apply to IC^^Qn, with the twisted W3 null-vector equations (143), which 
apply to @n^n- We will first do the comparison for general values of b, and then explain in more 
detail what happens in the particular limits 6 ^ 00 and b 0. 

The comparison for general b. To start with, the non-differential terms agree. This is actually 
a very non-trivial statement, as we started with complicated non-differential terms in eq. ( |4.35 ) 



an then generated more terms by twisting with 0„. The freedoms to choose the three parameters 
A, /u, u of Qn and to ignore terms belonging to ^2(2/) is a priori not sufficient to ensure the dozens 
of required cancellations, which nevertheless occur as can be seen in explicit calculations. These 
calculations use some helpful identities which are gathered in Appendix |A.1[ The existence of a 
simple twist which simplifies the differential equations obeyed by correlation functions involving 
many identical degenerate fields might well be a general phenomenon in conformal field theory, 
as we now see that it happens for the simplest degenerate field in theories with W3 symmetry, in 
addition to the already known cases of the two simplest degenerate fields in theories with Virasoro 
symmetry 

Let us then examine the term p-i^2 in eq- ( 143 ). Agreement with the corresponding term in 



eq. (4.22) would occur provided 



^■y^^lC-'^lclD2. (4.46) 

oy ^ y - Zi dzi 

It seems technically challenging to check this identity. But remember that our inability to explic- 
itly perform Sklyanin's change of variables for ^ does not contaminate the other terms in our 
equations, as we do know that the change of variables must be independent from the parameter 

b= (A;-3)-i 

Let us now examine the term p-Di. We would like this term to vanish modulo .^2(1/), as no 



such term is present in eq. ( 4.22 ). However, it is rather obvious that Di does not belong to T2{y), 
although it has quite a few remarkable properties. This is explained in detail in the Appendix |A.3| . 
As a result, the conjecture cannot hold for general values of b. 

The critical level limit b ^ 00. We notice that the term pi^i, which is responsible for the fail- 
ure of our conjecture, vanishes in the 6 — > 00 limit. Therefore, the conjecture has better chances to 
hold in that limit. To completely prove that it does, we still need to clear one subtlety with the term 
^D2. This term seems to vanish in the & — > 00 limit but actually it does not. This is because near 
b ^ 00 our correlation functions do not have finit limits. Rather, the Toda correlation function 

= {UaV-b-^w^{ya)U^V-bj,+b-^e^+e2)iz^)) inVOlvCS "hcavy" fields V_bj^_^b'He,+e2)iZi) 

whose momenta grow as b. On general grounds (see for instance [||]), it is therefore expected 
that ~ e''''^(^^'»r„ where S and Tn are 6-independent functions, and 5 depends only on 

{zi} and not on {ya}- The differential operator yiD2, which contains derivatives with respect to 
Zi, may yield a finite contribution when such derivatives act on e^^"^*^^^*^^. 

We should therefore check whether eq. ( [4.46 ) holds to the leading order in 6^ when acting 



on functions of the type ^'^^^^^^Tn- This is actually the case, because the only term in D2 with 
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Zi -derivatives is ^ • iF^I^' ^^'^ ^^^A;^ ^{{zi}) = ^5'({zi}). This completes the proof 
of the conjecture ( |1.4[ ) in the critical level limit b ^ oo 4^ k ^ 3. 

Notice that this b —>■ oo Umit is not sensitive to the twist function This is because the 



exponents A, /Li, 1/ (4.42) vanish in this limit so that — > 1. 

b— >oo 



The minisuperspace limit b ^ 0. In this limit, the discrepant term p-Di, which is responsible 



for the failure of the conjecture ( |1.4| ) for general b, grows larger. We may therefore obtain some 
insights on the reasons for this failure. 

As in the s£2 case, we will consider full correlation functions (with both holomorphic and 
antiholomorphic dependences) and use path-integral reasonings in 5^3 Toda theory. For full corre- 
lation functions, the conjecture reads 0„ ~ /C^|0„pO„. As in the si2 case, the transformation 
IC is 6-independent, r2„ is expected to have a finite limit, and the Toda correlation function 17„ 
behaves as f2„ ~ (rialli'^ V-b-^ioAya) 117=1 ^fe-i(ei+e2)(^*) / ^^ere Rn is 6-independent. 

b— >0 \ I 

Therefore $7„ simplifies in the 6^0 limit but, in contrast to the SI2 case, its leading behaviour 
does not reduce to a free field correlation function. This is because the simplified correlation func- 

r3n— 6 



tion ^Ha^i ^-b-iwi(ya) nr=i ^b-i(ei+e2)('2«)/ ^ocs not obcy momentum conservation, given 
the value 2Q = 2(6 + 6^^)(ei + 62) of the background charge in 5^3 Toda theory. However, mo- 
mentum conservation can be restored by inserting n — 2 screening operators l^-ig^. (See for 
similar reasonings and calculations in 5^3 Toda theory.) Thus, 

/3n-6 n n-2 \ 

,Zo ^ \ n ^-b-^u.Aya)ll^b-He^+e,){^^) H / ^'^^ ^b-^eA^i) ) ■ (4-47) 

^ \ a=l i=l i=l I 



This free correlation function is the product of the free correlation function ( [4.41| ), which we took 
as our ansatz for IG„|^, and an integral over xi, leading to 

\Qn?Cln Rn / \{cPxi J| \xi - ya\^ JJ '^^ " '^^l''^ 11 " ^^''""^ " '^"^"'^^^ 

^ £ a,t i,£ £^£' 

The integral in this formula is expected to be dominated by a saddle point, where the x^s are 
solutions of 

2 V + y — V = . (4.49) 

(Curiously, these are the Bethe equations for the SI2 Jaynes-Cummings-Gaudin model at infinite 
coupling and with spins it^ [23].) The dominant behaviour of the integral is expected to be of the 



2 2 



form 6^ as 6 ^ 0, with On a 6-independent quantity. This \6n\^ factor contradicts the existence 

of a finite Umit for as 6 — > 0, which follows from the conjecture. 

_ 1 

One may be tempted to modify the conjecture by adding a factor On '' to the twist function 0„. 
This would not only correct the leading behaviour in the b limit, but also make the conjecture 
compatible with global conformal symmetry. We have not mentioned global conformal symmetry 
until now because this subject is independent from the differential equations in terms of which the 
conjecture was formulated. It is however easy to see that the conjecture is incompatible with the 
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behaviour of correlation functions under scaling transformations {zi,ya 
the 6 — > oo limit. 

1 



{Xzi, Xya), except in 



However, adding the factor 6'„ would spoil the agreement between most terms of the KZ 
equations ( 4.22 ) and the twisted Ws null-vector equations ( [4.43 ), in particular the terms depend- 
ing on the spins jj. The modified conjecture would only hold at the level of the -independent 
dominant factors in the 6 — > limit, which would not be interesting. 



5. Conclusion 



The comparison of 5^3 KZ equations in Sklyanin variables (4.22) with Ws null- vector equations 



( |4.43[ ) does not support the conjecture ( |1.4[ ) in its general form. Nevertheless, the KZ equations 
are very similar to the null- vector equations: many terms agree nontrivially, and the disagreement 
is confined to a term which does not depend on the spins ji of the fields. This remarkable quasi- 
agreement makes it unlikely that a full agreement can be obtained by modifying the conjecture. 
In the critical level limit k —>■ 3 ^ b 00, the disagreement disappears and the conjecture 



(|1 .4|) is true. This limit plays an important role in the Langlands correspondence [24], which might 



possibly explain why the conjecture ( |1.4[ ) holds for s£2 and not for sis, and why in the si^ case 
it holds only in the critical level limit. Another hopeful source of insights is the recent work on 
conformal Toda theories [^], where the s£n>3 cases are understood to be qualitatively different 
from the 5^2 case. Of course, we already pointed out a significant qualitative difference, namely 
the failure of the s^s cubic field W'^{z) ( 2.1l| ) to obey the W3 algebra. It is not clear how this is 
related to our problem. 

Our results in the si^ case lead to natural conjectures in sij\[^3 cases, where we expect the KZ 
equations in Sklyanin variables to agree with Wn null-vector equations only in the critical level 
limit k —>■ N. Let us tentatively perform a counting of equations. There are ^N{N — 1) isospin 



variables on the Ihs of eq. (L4), and on the rhs we expect ^N{N — l)(n — 2) Sklyanin variables ya 



plus A^(A^ — 1) extra variables, which may be collectively included in the symbol U. Differential 
equations for the s£j\[ Toda correlation function which generalizes fi„ are obtained by eliminating 
i {N — 2) {N + l)n non-differential terms from the ^N{N — 1) (n — 2) null- vector equations. Thus, 
we have n — N{N — 1) differential equations. When it comes to /C0„fi„, we should presumably 
add an equation for each one of the extra variables, reaching n differential equations. This precisely 



the number of KZ equations for the Ihs ri„ of eq. (1.4). In addition, we have the same number of 
global Ward identities on both sides of eq. ( [L4] ), namely A^^ 



1. 
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A. A few technical results 



A.l Helpful identities 



The following identities are used in computing the non-differential terms of the operator = 
QnE2@n^ eq. (iA3). Some identites are written modulo terms in ^2(2/) ( |4.21 ), as indicated by 
the relation sign r^. All identities are proved by elementary manipulations, using observations of 



*e type ^ 



1 1 



tJa-Zi y~Zi 



{y-zxf 



1 

^ y-Vbivb- 



^y-Zi^ {y- ZjY 5Z r„ _ 



+ 



^ y-Vb{yb- Zif {y - Zif {y - zif y - Vb 



y-yh \,yb-Zi 



E 



1 1 



^y-yty- Zj yt - Zi 



^ [y - ^y-yby^y-Zi^ 

— V: ill — y-'lS /L^ n, _ I jL^ 1, _ 



E— 



1 1 



{y - ZiY" ^ y-yb 

1 y^ 

— r.A ^—^ 7/ — ' 



y- Zi 



(A.l) 
(A.2) 

(A.3) 

(A.4) 
(A.5) 



^ y-ybyb- Ziyb- zj {y - Zi){y - Zj) ^ y - yb 

1 1 1 ^ 1 1 

- yA'^ ' ^ 11 — 111 n — "y ■ * ^ iif — iiiX^ 

2 



y^ L_ 

u ^y- yb - z 



+ 



{y - ZiY (y - ZiY ^ y-yb y- Zi^ {y- ybY 



E 

b+c 



1 1 



1 



1 11 

+ 



y-ybyb-ycyb- Zi {y - Zif 2y- z, 



y -yb 



E^ 



1 1 



b^c 



y-ybyb-ycyc- Zi [y-zif 2y-z.. 



y ' 

^(y- ybY 

y ' 

b y^^^ 
1 1 1 1 1 

[y - ZiY ^ y -yb y - Zi ^ y - yb ^ yc- Zi 



2 1 1 

+ 



E^ 



y -yb 



(A-V) 

(A.8) 
(A.9) 



— ^ 1 1 ^ — ^ 1 ^ — ^ 1 ^ — ^ 1 

^ {y - ybY yb-yc ^ (y - ybY y-yc ^ (y - y^)^ ' 



E 

b^c 

E 

c^b 
dj^b 

E 

by^Cy^d 



b^c 

1 1 



1 



y-ybyb-ycyc - yd 



1 1 



1 



y-ybyb-ycyb- yd 



1 1 



1 



y -yby -ycy - yd 



1 1 



1 



1 1 ^ 1 

tf y-ybiyb- ycY 6 y - yby - ycy - yd 

b=f:c b^c^d 



y^ 1 1 1 y- 



1 1 



1 



y~yb (2/6 - ycY 3 y - yby - ycy - yd 



(A. 10) 



, (A.ll) 



(A. 12) 



E^l +E 

^y-yb I ^ 



(y - ybY 



3E^Er^3) 

^ y-yb'^ \y- ycY 
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A.2 A characterization of J^2{y) 



Here we will justify the characterisation (A.22) of the space ^2(2/) defined in eq. (121). 

For pedagogical reasons we will begin with the simpler problem of characterizing the space of 
permutation-symmetric functions of m variables {ya}- More precisely, given a function f{t, {ya}) 
which is permutation-symmetric in [ya], depends on an additional variable t, and is regular at 
t = ya, we want to determine whether /(y, {ya}) is actually permutation-symmetric although it 
apparently depends on y. This amounts to determining whether f{ya', {ya}) actually depends on 
the choice of a'. If it does not, then for any polynomial P{t) of degree m — 2 we have 

So we have transformed the m — 1 conditions /(yi) = f{y2) = • • • = f{ym) into the condition 
V„, S dt J^^9J^^\ = 0, which can then be evaluated by moving the integration contours, if the 
analytic properties of f{t) permit. 

Let us apply a similar reasoning to the characterization of J^2{y)- If f{y) G ^2{y), for instance 
/(y) = ( ^2 f{y) where fiy) is actually permutation-symmetric, then given any polynomial 
P{t) of degree n — 7 we have 



a' •'Va' 



dt Pit)^^^ —fit) = f(y) (p dt Pit)^^^^ — = . (A.15) 



Thus, to know whether f{y) G J^2{y), we only need to evaluate the left hand-side of this equality. 
To do this we can use the assumed analytic properties of f{t): namely, that it is meromorphic with 
singularities only at t = Zj, and goes to zero as t ^ 00. This implies 



which proves f{y) G ^2(2/) ^ {P, /) = as in eq. ( A.22| ). The reverse implication follows from 



a simple counting of variables: the space of polynomials of degree n — 7 has dimension n — 6, 
which is precisely the number of constraints which we expect for characterizing the space ^2(2/)- 

A.3 Study of the differential operator Di 



As explained in Section i.4, our conjecture (1.4) impUes the relation Di or equivalently 



Di G ^2{y), where Di is the first-order differential operator written explicitly in eq. ( 144 ). Here 
we provide a rigorous argument that this relation is not true, which implies that the conjecture 
cannot hold for general values of the parameter b. 

To start with, let us reduce the study of the first-order differential operator Di to the study of 
mere functions. The operator Di, like all our differential equations, is assumed to act on functions 
which are symmetric under permutations of the 3n — 6 variables {ya}- The space of such functions 
is algebraically generated by the 3n functions 

= ^log(ya-z,) , ai = ^-^ , Ti = Y,7-^-X2 ' i = l---n.{A.n) 

ya „ \ya ■^tj 



a a 
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Therefore, Di 

Dipi ~ 
Diai 

DlTi ' 



Dipi ~ DiCTj ~ DlTi ~ 0. Direct calculations show 



(y 



+ 



+ 



+ 



{y - ZiY {y - ZiY y 



6(7,; 



(y 



{y - z^Y 
1 



+ 



10 



(y 



z*)^ ^ y ■ 



+ 



4r, 



y 



+ 



{y - ziY tfAy- 



{y - Zi)^ 
4 

(y - Zi)'"" 



1 



y 



(A. 18) 
(A. 19) 

(A. 20) 

2 

|(A.21) 



So DicTj and Dir, do not manifestly vanish modulo J^2{y)- Let us however study them further. 
They may be considered as values at t = y of functions f{t) = f{t, {ya}, {zi}) which are invariant 
under permutations of {ya} but depend on the additional variable t. Let us consider the space of 
such functions, which we in addition assume to be meromorphic in t with no singularities besides 
t = Zi, and to go to zero as t ^ oo. Let us moreover introduce the space Vn-j of polynomials 



P{t) of degree n — 7. As we show in Appendix A.2, 

n „ 

f{y) € T2iy) 4^ VP G Vn-7, {P, /) ^ ^ ^ (it P{t) 



1=1 



Utiit 



Ut-i'it 



ya) 



-fit) = . (A.22) 



Then, explicit calculations yields 

.Uk^zi 



{P,Diai) = 27Ti- 



Zk) 



{P,DiTi) = 2Tri 



UkM^i 



-P'iz,. 



ya) 



Zk) 



Uai^i-ya) 



4P"izi) + 2(Ti + 6 



kj^i 



Zk 



P'{z.i) 



(A.23) 
(A.24) 



This explicitly demonstrates that Di ^ ^2{y)- 

However, Di still has remarkable properties with respect to the constant polynomial P = 1, 
namely (l,Di(Ti) = (l,DiTj) = 0. These non-trivial identities sensitively depend on the general 
structure of Di and on the particular values of A, ^, which determine its coefficients. This implies 
that, whereas arbitrary differential operators belong to T2{y) for n < 6, Di £ ^2(2/) for n < 7. 
The significance of these properties of Di is not clear. When combined with Dipi ~ 0, they 
suggest that Di ~ when applied to a special class of permutation-symmetric function of ya (and 
Zi), and one might wonder whether Qn^n actually belongs to this class. Given the freedom to 
choose y G {ya}, this would imply that Qn satisfies n — Q further differential equations. But Cl^ is 

not expected to satisfy any further differential equations besides the global Ward identities, whose 

7 

number is n-independent. So the supposition Di ■ Qn^n ~ certainly fails for n > 7, and so does 



our conjecture (1.4). 
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